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ABSTRACT 

We present a grid of nonequilibrium ionization models for the X-ray spectra from supernova rem- 
nants undergoing efficient diffusive shock acceleration. The calculation follows the hydrodynamics of 
the blast wave as well as the time-dependent ionization of the plasma behind the shock. The ion- 
ization state is passed to a plasma emissivity code to compute the thermal X-ray emission, which is 
combined with the emission from nonthermal synchrotron emission to produce a self-consistent model 
for the thermal and nonthermal emission from cosmic-ray dominated shocks. We show how plasma 
diagnostics such as the G'-ratio of He- like ions, defined as the ratio of the sum of the intercombination, 
forbidden, and satellite lines to the resonance line, can vary with acceleration efficiency, and discuss 
how the thermal X-ray emission, when the time-dependent ionization is not calculated sclf-consistently 
with the hydrodynamics, can differ from the thermal X-ray emission from models which do account 
for the hydrodynamics. Finally we compare the thermal X-ray emission from models which show 
moderate acceleration (~ 35%) to the thermal X-ray emission from test-particle models. 
Subject headings: cosmic rays - thermal emission: ISM - shock waves - supernova remnant - X-rays: 
ISM 



1. INTRODUCTION 

In young supernova remnant (SNR) shocks, the accel- 
eration of cosmic rays (CRs) leads to a softening of the 
equation of state in the shocked plasma. This comes 
about because the diffusive shock acceleration (DSA) 
process turns some non-relativistic particles into rela- 
tivistic ones and because the highest energy relativistic 
particles escape from the shock. Both of these effects 
lead to lower post-shock plasma t emperatures as well as 
higher post-shock densit ies (e.g., Uones fe Ellison! H99l 
IBerezhko~fc Ellison 1999|). The ionization of shocked gas 
at a particular time depends upon both the gas den- 
sity and the elec tron temperatur e , and it has been re- 
cently shown by iPatnaude et al.l (|2009ft that the time- 
dependent nonequilibrium ionization (NEI) and electron 
temperature, T e , are influenced by the DSA efficiency, 
cdsa- 

A number of young SNRs show some evidence for 
both nonthermal and thermal emi ssion behind the 
forward shock, in cluding SN 1006 dVink et all [200l 
iBamba et aH 12008ft . T vcho (jHwang et a l.l 12002ft. Kepler 
([Reynolds et al.ll2007t ). and recently Cas A (jArava et al.l 
12010ft . The thermal emission arises when the forward 
shock sweeps through the circumstellar medium (CSM) 
and h eats it to X-ray emit ting temperatures. As pointed 
out in lEllison et al.1 (|2007l ). the thermal X-ray emission is 
often considerably fainter than the nonthermal emission, 
but there exist examples where the thermal emission is 
as bright or brighter than any nonthermal e mission, such 
as in parts of RCW 86 (jVink et al.l [200l and N132D 
(|Xiao fc Chenll20M ). In contrast, in SNR RX J1713.7- 
3946, the lack of thermal X-ray emission is an important 
constraint on the ambient medium density, and we have 
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recently shown that if the TeV 7-ray emission is hadronic 
in origin, then copius thermal X-ray emiss ion would be 
obser ved with current X-ray observatories (jEllison et al.l 
12010ft . 

Here we extend our p revious work ([Patnaude et al.l 
12001 lEllison et all \20Wt) by coupling our models to a 
plasma emissivity code, producing thermal X-ray spec- 
tra and ionization timescale parameters (i.e. n e t) as a 
func tion of position and tim e in the SNR shock. Again, 
as 111 IPatnaude et~aTl (|2009ft . we limit our analysis to the 
region between the forward shock (FS) and contact dis- 
continuity (CD). In § 2 we present our new model and 
outline how we have coupled our hydrodynamical model 
to a plasma emissivity code. The inclusion of the hy- 
drodynamical evolution makes it difficult to compare the 
computed thermal emission to existing models, since the 
shocked, ionizing plasma is expanding adiabatically and 
the electron temperature is constantly evolving. To test 
the code, we present in § 2 a test where we follow the 
ionization in the case of a fixed density and temperature, 
which is equivalent to available NEI models in spectral 
fitting codes such as Xspec. 

The parameter space involved in producing the ther- 
mal and nonthermal emission as well as the hydrodynam- 
ical evolution is extensive and in § 3 we show a simplified 
set of models where we only vary a subset of parameters 
such as the shock acceleration efficiency. In this section, 
we show how the so called G-ratio, which is a measure of 
electron temperature, varies with acceleration efficiency. 
We find that since the increased acceleration efficiency 
generally leads to lower electron temperatures, increas- 
ing the acceleration efficiency results in higher G-ratios. 
In § 4, we compare our models to similar models where 
the X-ray emission is calculated not from the ionization 
state of each element, but from the ionization timescale 
which the shocked gas is expected to have at the end of 
the simulation. We find that the results differ, mainly 
at lower energies. Finally in § 4 we discuss the recent 
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result bv lBerezhko fc Volkl (120101 ) and show that the ab- 
sence of thermal X-ray emission does not naturally lead 
to the interpretation that the TeV gamma-ray emission 
is hadronic in origin. We argue that if the blastwave 
has only recently hit the circumstellar shell, then the 
shocked gas, though at X-ray emitting temperatures, will 
be underionized. W e show that our results and those of 
lEllison et al.l pOloh are consistent with the material be- 
ing at a lower density. 

2. EMISSIVITY CODE 

In previous papers (lEllison et al.ll2007t IPatnaude et al.l 
20091 lEllison et alJl2010l) 7 we outlined our model which 
coupled the NEI calculation to the cosmic ray accel- 
eration and hydrodynamics code. For each simulation 
cell and timestep, we calculate an ionization state vec- 
tor f(X l ), an electron temperature T e , and electron and 
ion densities, n e and n^. Combined with a chosen set of 
abund ances (here we choose those of lAnders fe Grevessel 
(|1989j )). these quantities are all that is required as in- 
put to produce the thermal X-ray emission. As in 
IPatnaude et al.l (|2009[ ). the ionization processes included 
in our NEI calculation include direct collisional ioniza- 
tion and collisional excitation followed by autoionization. 
Additionally, the recombination rates include radiative 
recombination and diclectronic recombination. We ne- 
glect Auger ionization as well as photoionization. In the 
current version of our code, we do not include ioniza- 
tion from super-thermal particles; in an io nizing plasma, 
these effects are expected to be small (jPorquet et al.1 
I2001afl . 

We have chosen to couple our models to an updated 
version of the Raymond-Smith emissivity code (RS93; 
IRavmond fc Smith! fl977t IBrickhouse et al.l[l995h . RS93 
differs from other available codes in terms of the number 
of included emission lines, a few line centroids, and indi- 
vidual line emissivities. These differences arc sufficiently 
numerous to make it difficult to compare our calculated 
spectra to other available models suc h as APED, but th e 
differences are well documented (c.f., iSmith et al.ll200i| ). 
Additionally, since our models not only follow the time 
dependent ionization, but also are dynamical in nature, 
it is not adequate to compare the calculated spectra to 
existing models such as nei and pshock in Xspec. 

To determine whether the modeled spectra are con- 
sistent with existing NEI calculations, as well as con- 
sistent with models in collisional ionization equilibrium, 
we performed the following test: we ran a simulation 
for a prescribed amount of time, and far upstream from 
the shock, we fix the electron temperature at a constant 
value of T e = 3x 10 6 K and keep the density fixed at n p $ 
= 3.0 cm~ 3 . That is, while in general the hydrodynam- 
ics computes the electron temperature and densities self- 
consistcntly behind the shock, for the purposes of this 
test we keep the electron temperature at a high, constant 
value far enough upstream of the shock so that we can 
evolve the NEI. This is equivalent to the single tempera- 
ture, single ionization age models found in Xspec, such 
as nei and vnei. We compare our computed spectrum 
to both what is computed in the Xspec model as well as 
a model for a plasma in collisional ionization equilibrium 
(e.g. raysmith or mekal). For simplicity, we limited our 
calculations to spectra that only include H, He, and Si. 
In Figure [TJ we show the computed and reference spectra 



for a shocked plasma that has been evolving for 500 yrs, 
as well as what the same spectrum looks like under con- 
ditions of collisional ionization equilibrium. In the top 
panel of Figure [U we also plot the evolution of Li- and 
He-like ions of silicon. 

At an age of 500 years, corresponding to an ioniza- 
tion age of ~ 5xl0 10 s cm~ 3 , at an assumed density of 
n p,o = 3.0 enr 3 , we observe emission lines from several 
atomic transitions. For instance, from the helium iso- 
sequence, we see line emission at 2.295 keV (Is 2 - ls4p), 
2.185 keV (Is 2 - ls3p) 1.866 keV (Is 2 - ls2p), and the 
multiplet transitions Is 2 - ls2p 3 P° at 1.85 keV and Is 2 - 
ls2s 3 S. Also, at lower energies we see numerous emission 
lines from the lithium to nitrogen iso-sequences, includ- 
ing multiplet transitions which lead to multiple spectral 
lines. 

In Figure El we plot the same model (T e = 3xl0 6 K 
and n Pi o = 3.0 cm~ 3 ) at an age of 2500 yrs. After 2500 
years, the same helium-like lines are seen as in the £snr 
= 500 year model, but the intensities have increased rela- 
tive to the lower ionization states. This is reflected in the 
top panel of Figure [2J where the lithium-like state of sil- 
icon is approaching the collisional ionization equilibrium 
value of ~ 0.1 (at T e 3xl0 6 K.) Likewise, in collisional 
ionization equilibrium at 3xl0 6 K, the He-like state is 
populated at ~ 0.87 (in fact, in collisional ionization equi- 
librium, Li- and He-like states are the dominant species 
at this temperature for silicon, with additional contribu- 
tions from the Be-like state; a negligible fraction of ions 
are in the H-like state), and the NEI conditions shown 
here are seen to be approaching it after 2500 years. 

In both test cases, the computed and reference NEI 
spectra are in good agreement, showing the same under- 
lying continuum and similar line emission, with differ- 
ences most noticeable in the emission lines between 200 
and 500 eV. We attribute these differences to small dif- 
ferences between the underlying atomic data used in the 
codes. The most obvious difference between the two is 
seen as the obvious lack of an emission line in the Xspec 
nei model at ~ 350 eV (Si XII n=2 to n=4 multiplet). 
In the tsNR = 500 year model, this line is quite bright, 
while in the higher aged model it appears to be weak- 
ening, relative to the surrounding lines (which are also 
weaker, relative to the continuum. 

The goal of this testing was to ensure that the NEI cal- 
culation in our model, in the limit of a single ionization 
age, single temperature nei model, is consistent with 
existing models. Realizing that there are differences in 
the underlying atomic data that will be manifest in the 
existence and intensity of individual emission lines (in 
particular, the nei model includes many more emission 
lines from He-like states, as seen in Figures [T] and [2]), 
our aim was to track the trending in the two NEI mod- 
els to ensure that they behave consistently against one 
another. As an additional test, we chose to compare our 
spectrum to that from a plasma in collisional ionization 
equilibrium. We find that, as expected, the computed 
NEI is in fact approaching the CIE spectrum expected 
for the temperature and density we are using for the test 
(jSmith fc Hughes Jl2010h . 



3. RESULTS 
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The cosmic ray hydrodynamics code contains an ex- 
tensive set of parameters which affect the SNR dynam- 
ics, emitted thermal and nonthermal broadband spec- 
trum, and relativistic particle populations. An exhaus- 
tive study of the effects that all these parameters have 
on the emitted spectrum is beyond the scope of this pa- 
per, but we note that combinations of parameters have 
been successfully used to produce emitted spectra from 
SNR RX Jf 7f 3.7-3946, ru ling out certain models for that 
SNR (jEllison et al.M2010Q . To simplify the process, we 
choose a set of models where we only vary the acceler- 
ation efficiency 6dsa and ambient medium density, and 
leave everything else fixed. We note that the applica- 
tion of these models to a specific SNR would require a 
very specific set of parameters, so we restrict ourselves 
to phenomenological studies for now. 

Our model assumes the following set of parameters: 
isNR = 1000 yr, E SN = fO 51 erg, and M e j = 1.4M Q . 
Additionally, we assume an exponential ejecta profile 
and an ambient magnetic field of I5/iG. Parameters 
which affect the ionization calculations and resultant 
thermal emission include the assumed abundances and 
the electron heating. W e assume cosmic abundances 
(j Anders fc Grevessell 19890 and heating via Coulomb col- 
lisions. The ambient medium has a temperature of 10 4 
K and is preionized at 10%. We vary the accelera- 
tion efficiency, edsa between 1 % (test particle) and 75%, 
and except where indicated otherwise, we chose ambient 
medium densities, n Pi o ~ 0.1, 0.3, 1.0, and 3.0 cm -3 . Our 
model also includes the ability to apply interstellar ab- 
sorption as well as combine the thermal and nonthermal 
emission. For clarity, we do not include any absorbing 
column, and except where noted, we neglect the contri- 
bution to the spectrum by the nonthermal continuum. 
We note, as in SNR RX Jf 7f 3.7-3946, that the underly- 
ing nonthermal continuum can be the dominant source 
of X-ray photons, so the models like those presented here 
are most appropriate for comparison with emission line 
fluxes for shocks that show a combination of thermal and 
nonthermal emission. 

The emitted thermal spectrum is a combination of 
brcmsstrahlung and two photon continua and line radi- 
ation. The formation of lines mainly occurs by electron 
impact of atoms or ions. The excitation can generally 
be broken down into excitation of outer-shell electrons, 
excit ation of inner-shell electrons, and resonant excita- 
tion ((Raymond k, Brickhouselll9960 . Additional contri- 
butions from radiative recombination to excited states, 
dielectronic recombination satellite lines and inner shell 
ionization are sometimes significant, especially in plas- 
mas out of ionization equilibrium. In the process of ra- 
diative recombination, the following occurs with 



+ e - Z +Z {n) + hv cont -> Z +Z (m) + hv mn , 

(1) 

where the ion Z + ( z+1 ^ captures an electron and is de- 
excited to the m-th excited state via emission of a pho- 
ton. In the process of dielectronic recombination, a free 
electron gains kinetic energy as it approaches an ion, then 
excites a bound electron to an energy level higher than 
the energy it had at infinity, so the free electron is cap- 
tured to form a doubly excited state of ion Z +z , 



£+(*+!) + e - (£+*)•* . (2) 

If auto-ionization occurs (the reverse of the capture pro- 
cess) the system returns to its original state and no re- 
combination takes place. Alternatively a fraction of the 
ions in the autoionizing state decays by spontaneous ra- 
diative transition of the inner excited electron to a state 
below the first ionization limit: 

{z +z y* (z +z y + hv, (3) 

where the stabilizing transition in the recombined ion 
Z +z results in the emission of a dielectronic satellite line 
of the parent transition in the rccombining ion 
Eventually, the singly excited state cascades down to the 
ground state with subsequent emission of photons. 

3.1. Emission Line Diagnostics 

To demonstrate a potentially useful diagnostic pre- 
dicted by the models that should be observable by Astro- 
H, we consider the n=2 — > n=l emission of He- like ions. 
This consists of a resonance (r) 2 1 P — > l x 5 line, a for- 
bidden (/) 2 3 S' — > l 1 S f line, and an intercombination (i) 
2 3 P -> l 1 ^ line. The intensity ratio (f+i)/r, called the 
G-ratio, varies with density, ultraviolet radiation inten- 
sity, temperature and ionization state. The shocked, 
swept-up material behind SNR shocks is considered to 
be at a low enough density (n e << 10 10 cm -3 ), and the 
UV radiation field is sufficiently negligible, that density 
dependent variations in the spectra are not observable. 
Therefore we will consider only the variation with tem- 
perature and ionization state. As a point of reference, the 
G-ratio decreases with increasing electron temperature 
l|Gabriel fc Jord an 1969). Finally, we note that there are 
a number of satellite lines (S) formed by transitions of the 
form ls2s 2 -> ls 2 2s, Is2s2p -> ls 2 2p and ls2p 2 -> ls 2 2p 
which lie at wavelengths near the i and f transitions. 
Since current and near term observations are unable to 
resolve their contributions to the spectrum, we define an 
observable ratio G' = (f+i+s)/r. 

The lines in the G' -ratio contain contributions from 
direct collisional excitation by electrons (/, i and r), ra- 
diative and dielectronic recombination of the H-like ion 
into the excited levels (/, i and r), dielectronic recom- 
bination of the He- like ion (s), innershell excitation of 
Li-like and lower ions (s) and innershell ionization of Be- 
like and lower ions (s). 

Since the resultant emitted thermal X-ray spectrum 
can be quite complicated, containing emission lines from 
H- and He-like states as well as a forest of emission lines 
from intermediate charge states (e.g. as was demon- 
strated in Figures [T] and [2]) , we focus on the bulk char- 
acteristics of the emitted spectrum. To do this, we sum 
the emitted spectrum from the contact discontinuity to 
forward shock and investigate the global characteristics 
of the thermal emission as a function of both accelera- 
tion efficiency, and to a lesser extent the ambient medium 
density. 

In Figures [3]- [6] we plot the G'-ratio for the He-like 
states of oxygen, neon, magnesium, and silicon. Because 
excitation cross sections for the forbidden and intercom- 
bination lines decrease with energy, relative to that of 
the resonance line, and because the cross sections for the 
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forbidden and intercombination lines have strong reso- 
nance contributions near the excitation threshold, G' can 
be thought of as a diagnostic of the electron tempera- 
ture. In addition, the dielectronic recombination satel- 
lites are excited by electrons well below the threshold for 
exciting the r line, but since they are observed near the 
triplet state (where they remain unresolved with current 
instrumentations), their contribution to the flux should 
be included, thus weakening the overall dependence upon 
temperature. The G'-ratio also depends upon the ion- 
ization state of the gas, since the triplets arc populated 
by recombination from the H-like states, and innershell 
excitation and ionization of lower charge states also con- 
tribute. The interplay of temperature and ionization 
state requires a fully time-dependent ionization model 
such as we present here, and this ratio will be an im- 
portant diagnostic tool for upcoming missions such as 
Astro- H. We note that we have so far neglected exci- 
tation and ionization by the nonthermal electrons (e.g. 
iPorquet et al.1 (|2001aD ). and that some improvements to 
the a tomic rates in o ur code remain to be implemented 
(e.g., IPorquet et al.l (|2001bD ). but that we expect the 
tren ds of the predict i ons to be reliable. 

In IPatnaude et al.l (|2009D we showed that the electron 
temperature varied roughly inversely with the accelera- 
tion efficiency, that is, as acceleration efficiency increases, 
the electron temperature decreases. Therefore, in order 
to understand how the G'-ratio changes, we plot it versus 
the acceleration efficiency, rather than electron temper- 
ature. In Figure [3J we show the variations in G' of oxy- 
gen with acceleration efficiency, at constant density. At 
high ambient densities G' behaves as expected. That is, 
G' increases with increasing efficiency (or, equivalently, 
increases with decreasing temperature). However, the 
lowest density ambient medium model shows the oppo- 
site trend, except at the highest acceleration efficiencies, 
where G' is seen to rise rapidly. Similarly, in neon, as seen 
in Figure 0] the G'-ratio decreases with increasing effi- 
ciency for n am i, = 0.1 cm" 3 . However, at a density of 0.3 
cm" 3 , the trend with acceleration efficiency is consistent 
with the high density results for oxygen. A similar result 
is seen for magnesium, while in silicon, the variation in 
G' at n am b = 0.1 and 0.3 cm -3 is opposite to that seen in 
the higher density models (compare the upper two panels 
of Figure [6] to the bottom two). The decline in G' with 
acceleration efficiency in the lower density models can 
be understood as an effect of the underionization of the 
plasma. An increase in acceleration efficiency leads to a 
higher compression ratio, and since the ionization rates 
of ions below He-like are not very sensitive to tempera- 
ture for these ions at the temperatures of these models, 
the higher acceleration efficiency models reach higher He- 
like ionization fractions. This increases the strength of 
the resonance line relative to the inner shell excitation 
lines. 

3.2. Global Properties 

To better understand the variations seen in G' with 
cdsa, we plot in Figures [7H5] the variation in G' as a 
function of weighted electron temperature, at constant 
ambient medium density. The weighted electron temper- 
ature is over the entire region between the forward shock 
and contact discontinuity, and we weight the electron 
temperature in each hydrodynamical cell by the emission 



measure of the cell (i.e. <T e > = EM" 1 E^TJ x EM 1 )). 
The emission measure may not be the best quantity to 
weight the temperature by, and it may be more mean- 
ingful to weight the temperature by the power in the 
He-like triplet, but that will result in a weighted tem- 
perature which varies with element, such that the aver- 
age temperature as derived from the He-like oxygen will 
be different from the average temperature as derived by 
silicon. Using the emission measure to weight the tem- 
perature allows us to compare how G' varies with T e for 
each element, and also compare how it varies between 
elements. We note that the emission measure in each 
cell is the only quantity related to the emitted spectrum 
which will not vary with element . 

In Figure [7] we plot G' versus T e for O, Ne, Mg, and 
Si, at constant ambient medium density. The emission 
measure weighted average electron temperature for the 
shocked CSM varies from - 7xl0 6 K (~ 600eV) in the 
highly efficient models to > I0 7 K (~ 850 eV) in the 
test particle models. The curves are complicated; in 
the top panel, O is seen to decrease initially with in- 
creasing < T e > before increasing above 9xl0 6 K. How- 
ever, everywhere the resonance line intensity is well above 
the sum of the forbidden, intercombination, and satellite 
lines (G' < 1.). In neon, again the resonance line is 
brighter than the other contributions, while in magne- 
sium and silicon, the flux from the forbidden, intercom- 
bination, and satellite lines is greater than for the res- 
onance line, and increases with increasing temperature 
(i.e. decreasing acceleration efficiency). 

Generally, the G'-ratio is sensitive to temperature; at 
high temperatures, all states are collisionally excited, 
but at lower temperatures, dielectronic recombination 
becomes important. The triplet states will be preferen- 
tially populated through recombinations simply because 
of their greater statistical weight, thus increasing the in- 
tercombination and forbidden line strengths relative to 
the resonance line. This is seen in the low temperature 
(high acceleration efficiency) end of Figure where the 
G'-ratio is highest. This is not, however, the case for the 
lower density model. Here, the G'-ratio is generally seen 
to be highest in the test particle model. 

To understand the differences between Figure [7] and 
we plot in Figure [9] the average charge state for silicon 
and oxygen. Two trends are immediately clear from this 
plot. First, in the low density case, the average charge 
state for silicon is well below the He-like state, but the av- 
erage charge state increases with acceleration efficiency. 
Additionally, a higher ambient medium density results in 
a higher average charge. This is because the higher den- 
sity behind the shock increases the collisional ionization 
rate, and this dominates over the temperature depen- 
dence of the ionization rates of Si XII and lower ions. At 
the highest densities, the average charge is much higher, 
and it decreases with increasing efficiency. In this case, 
the decline in temperature with acceleration efficiency 
decreases the ionization rate of Si XIII, and this dom- 
inates over the increase in density because of the high 
ionization potential of He-like Si XIII. As seen in Fig- 
ure [9] similar trends are seen in oxygen. Here, however, 
even at the highest efficiencies in the high density cases, 
oxygen is almost fully ionized. 

Equating these results to Figures [7] and means that 
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in the low density case (Fig. [7]), the lower density im- 
plies less collisional ionization and less efficient Coulomb 
heating, so that the plasma stays less ionized. More pre- 
cisely, in the low density models, the He-like states are 
not reached, and forbidden line emission is dominated 
by inner shell excitation of lower ions. In higher density 
models, the He-like states are reached, and the expected 
dependence with temperature is seen (c.f. Fig. [8]). We 
stress that in the low density models, some elements, 
such as silicon, will be highly underionized and the He- 
like triplets will be quite faint. These lines will be hard to 
detect, and when combined with any nonthermal contin- 
uum emission, they might be undetectable with current 
instrumentation. 

We end this section with Figure [lOl a plot of the vari- 
ation in G', average charge state, and electron temper- 
ature as a function of ambient medium density, at con- 
stant acceleration efficiencies. This plot shows how these 
values vary between the test particle case and a 50% 
efficient shock. As expected, the average electron tem- 
perature in the test particle case is higher than in the 
efficient case, but the average charge states are very sim- 
ilar. The biggest differences come from the ratio of G", 
where in the low density models, the G' -ratio is higher 
in the test particle case than in the efficient case, while 
above a density of < 1.0 cm -3 , the relation is reversed, 
with the G' -ratio in the efficient case is higher. From this 
plot, it is clear that at identical ambient medium densi- 
ties, similar charge states can result between the efficient 
and inefficient models, while the temperatures can differ 
by ~ 30% and the line ratios can differ by ~ 10%. For in- 
stance, for an ambient medium density of 1.0 cm , the 
weighted temperature in the 50% efficient model is 30% 
lower than in the test particle case. However, the ioniza- 
tion state for silicon here is virtually identical, but the 
G 1 -ratio in the efficient model is 10% higher than in the 
inefficient model. At temperatures which can differ by as 
much as they do in Figure [TU1 what we see here is that at 
constant ambient medium density but differing efficien- 
cies, different lines are seen at different intensitcs. Since 
the line ratios and underlying electron temperatures are 
different between the efficient and inefficient models, the 
shape of the resultant spectra will differ. 

4. DISCUSSION 

The models presented here differ from existing avail- 
able models that compute the X-ray emission from 
shocks in several ways, but the relevant difference here 
is that the nonequilibrium ionization is followed simulta- 
neously with the shock dynamics and particle accelera- 
tion. The ionization vector is then passed directly to an 
cmissivity code to compute the thermal X-ray emission. 
Since the NEI is evolved simultaneously with the hydro- 
dynamic evolution, our calculation includes changes to 
tfie NEI as a result of adiabatic losses, or as is the case 
in an efficient model, increased compression behind the 
shock resulting from efficient shock acceleration. 

To illustrate this effect, in Figure [11] we compare the 
thermal X-ray emission from models where we follow 
the time-dependent ionization with the hydrodynamics 
(black curves) to the thermal X-ray emission from mod- 
els where we calculate the ionization structure from the 
final ionization age (r) and electron temperature (red 
curves) . Here, the final ionization for each cell age is de- 



fined as the r at the end of the simulation. Then, for 
comparison, we re-run the NEI calculation for the final 
electron temperature and the time equivalent to the age 
of the SNR at the end of the hydrodynamical simula- 
tion. For these models, we do not include any efficient 
particle acceleration, as that will only complicate the in- 
terpretation of these results. Computing the ionization 
structure in each grid cell from the final r is similar to as- 
suming a single r, single T e model for each cell in Xspec 
(e.g. the nei model), except that in our computation 
of r and T e in the hydrodynamics, we still account for 
adiabatic cooling in the shocked gas, since these values 
were computed initiall y from the hydrody namics. This 
approach is similar to lEllison et al.l (|2007l ). As seen in 
Figure QTJ there are differences between the two resultant 
spectra. In the low density (n p .o = 0.1 cm~ 3 ) model, the 
calculation of the emitted spectrum from the final ion- 
ization age slightly ovcrpredicts the final thermal X-ray 
emission, most evidently seen around the He-like lines of 
magnesium and silicon, as well as from the emission of 
L-shell iron. Not surprisingly, the shape of the under- 
lying continua are consistent. This is expected because 
the shape of the continuum is determined by the elec- 
tron temperature, which are identical between the two 
models. 

In the higher density (n Pi o = 1-0 cm" 3 ) model, the dif- 
ferences are much larger, particularly at low energy. In 
particular, the ionization state and thermal X-ray emis- 
sion computed from the final ionization age and tem- 
perature undcrprcdict the Fe-L emission. We interpret 
this difference as being a direct result of not following 
the ionization explicitly; that is, when computing the 
ionization state from a single density and At, interme- 
diary ionization states are not correctly populated, since 
the calculation assumes that the shocked gas has gone 
from a cold (10 4 K) gas to hot (several 10 6-7 K) in 
a single timcstep. We note that were a spectrum like 
this fit with existing models, it would appear that the 
shocked plasma is overabundant in metals such as argon 
and calcium. Figure Qj] clearly demonstrates the need to 
compute the ionization self-consistently with the shock 
hydrodynamics (regardless of whether one considers the 
effects of efficient diffusive shock acceleration). 

4.1. Efficient vs. Inefficient Models 

In Figure [l2l we plot the observed differences between 
the thermal emission from a test-particle model and the 
combined thermal and nonthermal emission from an ef- 
ficient model. The black crosses in Figure [T2l are sim- 
ulated data for a 50 ksec Chandra ACIS-S observation 
of an SNR with £dsa = and n p fl = 0.3 cm -3 . The 
blue curve is the nonthermal emission, the red curve is 
the thermal emission and the black curve is their sum. 
In the efficient model, we assume £dsa = 35% and n Pi o 
= 0.3 cm~ 3 . As expected, at the high energy end, the 
nonthermal emission dominates the shape of the spec- 
trum. However, below about 2 keV, where the the ther- 
mal emission is comparable to the nonthermal emission, 
differences are still seen. 

Below 1 keV, the thermal emission is ~ an order of 
magnitude above the nonthermal emission. It is in this 
regime where the differences between the thermal emis- 
sion from an efficient model and a test partcle model are 
most apparent. In particular, the emission from oxygen 
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is seen to differ between these two models, and is reflected 
in the fit residual, shown in the lower panel, which shows 
that the modeled thermal emission from the efficient case 
does not accurately describe the thermal emission in the 
test particle case. This suggests that the emission lines 
contribute at differing levels between the efficient and 
inefficient models. 

4.2. Comparisons to Recent Calculations for Limits on 
the Thermal X-ray Emission in RX 31113.1-39^6 

Recently. IBerezhko fe Volld (|2010l ) used models from 
lHamilton et al.l (1 1 9831 ) to constrain the origin of the TeV 
7-ray emission in RX J1713. 7-3946. They argued that, 
based on fits to the broadband nonthermal continuum, 
the TeV 7-ray emission is hadronic in origin, with ~ 35% 
of the shock kinetic ene rgy being dep o sited into nuclear 
cosmic rays . However, lEllison et al.l (|2010( ) argue that 
if the TeV 7-ray emission were hadronic in origin, ther- 
mal X-ray emission would be detected, though the model 
used there is for a constant ambient medium density of 
0.2 cm -3 , for the hadronic scenari o. In the model de- 
veloped bv IBerezhko fe Volld (|2010D . the SNR blastwave 
has been expanding through a low density (rt# w 0.0008 
cm" ) bubble and has recently encountered the bubble 
wall, with n h = 0.25 cm -3 at present. From their model, 
we estimate that the shock began encountering the bub- 
ble wall ~ 800 years ago, as determine d by the curve for 
N„ in the top panel of their Figure 1 (|Berezhko fe Volld 
2010). They assume that the thermal X-ray emission 
as det ermined for the Sedov solution bv lHamilton et ahl 
(|1983f ) is reduced by an amount equivalent to the ratio of 
the bubble emission meas ure to the Sedov emissio n mea- 
sure, here equal to 0.46 (jBerezhko fe Volld 120101 ). Ad- 
ditionally, t hey assume a value o f 77 = 8xl0 49 , where 
T) = n 2 EsN (jHamilton et al.lll983l) . Their chosen value 
of 77 corresponds to an explosion energy of 1.3 xlO 51 erg 
and an ambient medium density of 0.25 cm" 3 . They use 
the value of 77 and an assumed electron temperature of 
1 kcV to estimate the thermal X-ray emission from the 
recently shocked bubble wall, and they find that at 1 
kcV, the thermal X- ray emission is approxi mately half 
the nonthermal flux (jBerezhko fe Volkll2010!) . 

Since the blastwave hit the bubble wall ~ 800 years ago 
and the electron temperature is ~ 1 keV, the shocked 
plasma will be underionized. This is probably best il- 
lustrated in Figure I10[ where for an ambient medium 
density of 0.25 cm -3 , the average charge state of sili- 
con is quite low, and in fact since the model shown in 
Figure [TO] is for a constant density model, it represents 
an upper limit on what the ionization state would be in 
IBerezhko fe Vo lk (2010). Thus, the shocked plasma from 
the bubble wall will likely be underio nized. 

Basically, IBerezhko fc Volld (|2010t) use a model for 
thermal X-ray emission from a Sedov-like SNR (of age 
1600 yr and distance 1 kpc, but renormalized to match 
the emission measure of a shocked shell) to the emission 
from a recently shocked thin shell undergoing efficient 
particle acceleration. As we showed above, however, The 
emitted spectra in the two scenarios are expected to be 
markedly different. We stress that in order to constrain 
the origin of the TeV gamma-ray emission, whether it 
arises from either leptonic or hadronic sources, the time- 
dependent ionization and thermal X-ray emission must 



be calculated in a self-consistent manner with the shock 
acceleratio n and hydrodynami cs, as was done for a sim- 
ple case in lEllison et al.l ((20T0I) . 

As was shown in previous sections, the increased com- 
pression from the shock acceleration increases the line 
intensities (from oxyge n and neon, for instance), so 
IBerezhko fc Vo lk (2010j) will underestimate the contri- 
butions o f line emission to th e total flux. This is consis- 
tent with lEllison et ahl (|2010l) where both a low ambient 
medium density and a high clcctron-to-proton ratio at 
the maximum energies were required in order to hide 
any emission lines and still fit the shape of the GeV-TeV 
continuum. 

Since the shock acceleration does has an effect on the 
ionization balance and thus the emitted thermal X-ray 
spectrum, using thermal X-ray emission models which do 
not account for increased ionization or compression when 
trying to understand the origin of GeV-TeV emission is 
not accurate. This will be of particular importance when 
the blastwave is moving through complex environments 
such as circumstellar shells, where the density can change 
by a few orders of magnitude over relatively small spatial 
scales. In SNR RX J1713.7-3946, the absence of ther- 
mal X-ray emission rules out a hadronic origin to the 
7-ray emission in the case of expansion in to a uniform 
circumstellar medium (jEllison et al.|[20lol ). Additional 
modeling will be required to assess whether or not the 
hadronic origin to the TeV emission is a viable picture 
under other, more complex conditions. 

5. CONCLUSIONS 

We have presented an extensive grid of models where 
we study the effects on the emitted thermal X-ray spec- 
trum behind SNR shocks. To accomplish this, we have 
coupled our CR-Hydro+NEI calculation to an updated 
version of the Raymond-Smith plasma emissivity code. 
In studying the global properties of the shocked plasma, 
we find that in the low density limit, the plasma temper- 
ature diagnostic, the G' -ratio, decreases with decreas- 
ing temperature (increasing acceleration efficiency) for 
higher Z elements such as magnesium and silicon. This 
is opposite to the behaviour at higher densities, where 
the G'-ratio increases with decreasing temperature, as 
expected. We find that the reason for this behavior is 
that in the lower density models, the electron tempera- 
tures are lower, but more importantly, the ions do not 
reach the higher ionization states. Additionally, we find 
that in low density models, the average charge states 
increase with increasing efficiency, while at higher densi- 
ties, the charge state is seen to decrease. This is because 
in the low density model, the higher efficiency leads to 
higher post shock densities and therefore more collisional 
excitation, while in the high density models, the shocked 
plasma does not reach as high a temperature in the effi- 
ciently accelerating models, so the ions are less ionized, 
even though the post shock density is also high. The 
differences in average charge state between the low and 
high density models are thus viewed as a result of the 
temperatures in the low density model generally being 
higher than in the high density model. 

Based on our results, we conclude that it is not suffi- 
cient to calculate the ionization state of the plasma as 
a "post-processing" step. While in lower density models 
the differences arc small, at higher densities, the calcu- 
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lation will undcrprcdict the emission from intermediate 
charge states, particularly emission from Fc-L. We specu- 
late that these differences might be even larger in metal 
rich ejecta where the number of free electrons will be 
large. Finally, we note that the shape of the emitted 
spectrum in an efficient model, even when folded through 
a CCD-resolution response, is fundamentally different 
than the emitted spectrum from a test-particle model. 
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Fig. 1. — Top Panel: Evolution of ionization fraction for Si +11 and Si +12 in a single temperature noncquilibrium ionization calculation up 
to tSNR = 500 yr. The solid and dashed red lines correspond to the collisional ionization equilibrium values for Si+ n and Si +12 at T e = 3 X 
10 6 K. Bottom Panel: Thermal X-ray emission line spectrum for a single temperature NEI model at an ionization age of 4.7x 10 s cm , 
which corresponds to tsNR = 500 yr and n Pi o = 3.0 cm -3 . We only include emission from H, He, and Si. The black curve corresponds to 
our calculation and is compared directly to an equivalent Xspec nei model at the same ionization age and temperature, which is shown as 
the red curve. The blue, curve corresponds to the emission from a plasma at the same temperature in collisional ionization equilibrium. 
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Fig. 2. — Top Panel: Evolution of ionization fraction for Si +11 and Si +12 in a single temperature noncquilibrium ionization calculation 
up to tsNR = 2500 yr. The solid and dashed red lines correspond to the collisional ionization equilibrium values for Si+ n and Si"*" 12 at T e 
= 3 X 10 6 K. Bottom Panel: Thermal X -ray emission line spectrum for a single temperature NEI model at an ionization age of 2.4x10 
s cm~ 3 , which corresponds to tsNR = 2500 yr and n p> o = 3.0 cm~ 3 . We only include emission from H, He, and Si. The black curve 
corresponds to our calculation and is compared directly to an equivalent Xspec nei model at the same ionization age and temperature, 
shown here as a red curve. The blue, curve corresponds to the emission from a plasma at the same temperature in collisional ionization 
equilibrium. 
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Fig. 4. — G'— ratio for He-like Neon in models with n p fi = 0.1, 0.3, 1.0, and 3.0 cm 3 versus acceleration efficiency. The G-ratio is 
integrated over the entire region between the forward shock and contact discontinuity. 




Fig. 5. — G'— ratio for He-like Magnesium in models with n Pi o = 0.1, 0.3, 1.0, and 3.0 cm 3 versus acceleration efficiency. The G— ratio 
is integrated over the entire region between the forward shock and contact discontinuity. 




Fig. 6. — G'— ratio for He-like Silicon in models with n Pi o = 0.1, 0.3, 1.0, and 3.0 cm 3 versus acceleration efficiency. The G-ratio is 
integrated over the entire region between the forward shock and contact discontinuity. 
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Fig. 7. — G'— ratio as a function of emission measure weighted electron temperature for oxygen, neon, magnesium, and silicon, for an 
ambient medium density of 0.1 cm -3 . 
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Fig. 8. — G'— ratio as a function of emission measure weighted electron temperature for oxygen, neon, magnesium, and silicon, for an 
ambient medium density of 1.0 cm -3 . 




Fig. 9. — Average integrated charge states for silicon (left hand panels) and oxygen (right hand panels) as a function of density and 
acceleration efficiency. The higher density models show some scatter but, like the lower density models, there appears to be a clear trend 
with acceleration efficiency. 
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Fig. 10. — Top: The ratio of G' for He-like silion in a 50% efficient model versus a test-particle model for a range of densities between 0.1 
— 2.0 cm -3 . Middle: The average charge state behind the shock for the test-particle and efficient model. Bottom: The emission measure 
weighted electron temperature. 
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Fig. 11. — Top: Simulated Chandra ACIS-S X-ray spectrum from a model with n Pi o =0.1 cm 3 where the ionization is calculated sclf- 
consistcntly with the hydrodynamics (black curve) versus the ionization being calculated after the simulation, based on the final ionization 
age and temperature (red curve). Bottom: Same as in the top panel, but for an ambient medium density of n p .o = 1-0 cm -3 . In both 
cases, we assume cdsa = 1.0%. 
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Fig. 12. — Comparison between the emitted thermal X-ray spectrum from a test-particle model versus the combined thermal and 
nonthermal spectrum from a 35% efficient model, with n Pi o = 0.3 cm -3 . The data correspond to the simulated test-particle spectrum, the 
black curve corresponds to the combined thermal and nonthermal spectrum from the efficient model, the blue curve corresponds to the 
synchrotron component contribution to the model, and the red curve corresponds to the thermal contribution. The bottom panel shows 
the difference between the thermal emission from the test particle model (the data) and the model (here, the thermal emission from the 
efficient model as well as the contribution from the continuum.) 



